Slack Dynamics on an Unfurling String 



J. A. HaniiE^ and C. D. Santangekfj] 

Department of Physics, University of Massachusetts, Amherst, MA 01003 
(Dated: February 6, 2012) 



We consider the phenomenon of a noisy arch growing on a rapidly deployed thin object contacting 
a rigid plane. We present a qualitative model for the growing structure involving the amplification, 
rectification, and advection of slack in the presence of a quasi-steady stress field, validate our assump- 
tions with numerical experiments, and pose new questions about the spatially developing motions 
of thin objects. 
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The deployment of a thin object from a roll or pile is a canonical event, important for textile manufacture P~H3], 
tethered satellite control [U-[6], and successful use of anchor rodes, drogues, watch chains, and garden hoses. Despite 
this ubiquity, we were recently surprised by an observation of these dynamics, as embodied by a rapidly straightened 
chain on a tabletop |7 . The surprise in question is displayed in Figure [T One free end of an orderly monolayer of 
chain is rapidly pulled along a table, leading to the formation of a slow y growing, wobbling, noisy arch near the 
pick-up point. The phenomenon is not unique to chains, and may be observed in strings or ropes [5] and other 
quasi-one-dimensional objects. We will use the words "chain" and "string" interchangeably in what follows. This 
letter details some results that arise from consideration of the simplest possible system resembling that of Figure [T] 
which we believe shed some light on the phenomenon of the arch. 



FIG. 1. Still images from a video [7] showing configurations of a chain tightly woven on a table before (top) and during (bottom) 
a straightening process. A vertically-oriented arch is apparent in the latter image. Chain links are ~2 mm, the amplitude of 
the initial condition is ~5 cm, the imposed velocity of the pulled end is ~8 m/s. 

Given a time-dependent curve X(s,t) parametrized by arc length s, the wave equation 

M a t 2 x = d s {ad s x) (f ) 

describes a balance of inertia and line tension in a string of uniform mass density fi. The stress a(s,t), morally 
equivalent to the pressure in a one-dimensional fluid, is a multiplier field enforcing the metrical constraint c^X-c^X = f 
[51118]. This equation is in principle derivable as the zero-radius limit of an elastic rod, the continuum limit of a chain 
of rigid links or, more abstractly, from a Lagrangian living on a curve. Taking an arc length derivative of ([I]) and 
projecting along the unit tangent 9 S X yields the auxiliary equation 

2 s a - ok 2 = -fid t d s X ■ d t d s Tt , (2) 

the solution of which serves to impose the constraint [9l IT3HT8] . The "screening potential" of this Poisson-like equation 
is the squared curvature k 2 = d 2 H. • c^X, while the "source term" on the right hand side has been rewritten in terms 
of the magnitude of the rate of change of the unit tangent using c^c^X • 9 S X = 0. From this form, it appears that the 
stresses induced by inertial motions of the string will be tensile unless one imposes unfavorable boundary or initial 
conditions. Positivity of the stress is proven for one fixed and one free end [18] and the proof may be extended [19j . 
Hence, we will not generally be concerned with unstable evanescent dynamics that might arise from imaginary wave 
speeds in ([!]). 
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Relevant to our problem are a class of solutions to ([!]) and ([2| for which the curve motion is purely tangential, 
9 t X = Td s ~K. Here, X is arbitrary and T and a — [iT 2 are uniform along the curve [20]. The shape X itself is simply 
one of two solutions, namely that which travels at speed — T with respect to the Lagrangian frame of ([!]) moving at 
speed T along the curve. It is thus stationary in the laboratory frame; this and the linear dispersion inherent when 
a is uniform are used to great effect in certain kinetic sculptures [2T] [52] . 

The addition of a uniform gravitational term /ig to the right hand side of ([!]) breaks the degeneracy in X and 
may generate compressive stresses. The remaining stationary-shape solutions of this augmented equation are the 
catenaries |17l 123) with the tangential velocity T determining the tension, not the shape. Given that compression 
(<7 < 0) corresponds to catastrophic failure, we can naively assume inertial stabilization of a symmetric catenary arch 
a a negative constant, bearing a tension a = fj, [T 2 + ||g|| Vo 3 « 3 _±i \ as j on g as ^he length parameters of 



~~ a 2 s 2 + l ' 

the shape fall within a circle whose radius is the natural length scale defined by the velocity and gravity: (^) 2 + s 2 < 

(li^ll) ' ^ ^ s P om ^ ^ * s wor th recalling the scales involved in Figure |lj The terminal "free stream" tangential 
velocity T ~ 8 m/s, g is Earth gravity (~ 10 m/s 2 ), and k > ~ 10 /m for the duration of the experiment. Hence, 
the ratio of inertial to gravitational accelerations ^# is at least on the order of fifty to a hundred. So gravity is 
neglectable here, due to the large curvature that spontaneously emerges early in the process. The natural time scale 
of about one second (jrar) is equally irrelevant. Adding a rigid motion with comparable velocity to the analysis — a 
wobbling perpendicular to the catenary's plane, for example — does not change the order of magnitude of the relevant 
terms. However, given sufficient growth time, or a stiffer curve with a bending or twisting length scale that could set 
a lower bound on this initial emergent curvature, gravity may come into play. We will also neglect drag on the chain 
by air and table. 

Returning to our original equations ([I]) and ([2]), let's consider an initially planar curve X with a small perturbation 
<56b oriented perpendicular to the osculating planes of X, in this case simply the plane containing X. To first order 
in the perturbation, we have the decoupled equations 

Hd 2 db = d s (ad s Sb) , (3) 
d 2 Sa - ban 2 = , (4) 

such that 8a = 0, and the height function Sb moves through a stress background determined entirely by the planar 
dynamics. Modulo the periodic component of the motion induced by the initial chain layout, material particles take 
heteroclinic trajectories between two arbitrary stationary shapes with tangential velocities of zero in the pile and the 
imposed pulling velocity downstream. Although we might expect velocity and stress discontinuities at the moving 
front, we may also clearly observe a smooth transition region for the velocity [7], which presumably corresponds to 
a smooth transition in stress. The dynamics as a whole appear to be nearly steady; even the arch grows slowly and 
remains within a nearly constant distance from the front. Thus, with an imposed pulling velocity T, we introduce the 
traveling wave variable r\ = s + Tt, presume some a(rj), and approximate the slowly evolving dynamics Q with the 
equation 
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dt sb = a M - jg dSb 

2/zT n 



(5) 



This is a short-time description of b(r), t) that presumes smallness of second order time derivatives, in a frame Eulerian 
with respect to the curve but not the laboratory. It admits general solutions, for which the characteristics converge 

at infinite time, of the form Sb = f (t + J n drj' a ^yL^ T i ) , f an arbitrary function. Figure J^J compares the evolution 

of a wave packet under ^ and the integrated form of ^ for a tanh-shaped a(rj) that interpolates between zero 
and the free stream stress fiT 2 . Agreement is good for early times, when the packet is "squeezed" behind the free 
stream. The disturbance eventually reaches a nearly stationary configuration via a complicated modulated motion, 
not captured by the first order transport equation §5§ . Note that the area under the height function does not represent 
any "energy-like quantity to be conserved" [17] in this transport equation, so that in this approximation, squeezing 
the packet does not increase its amplitude. 

The initial squeezing is best expressed in a new variable, the "slack" I = y (d s 5b) 2 taken up into the height function. 
In terms of I, the transient dynamics §5§ take the form 

V(r?)-/iT 2 



^=H^V^J- (6) 

Given a monotonically increasing a(rj), Q is comprised of an amplification and an advection term. These are 
the ingredients for a convectively unstable system 2 1 27 . However, in contrast to the open flow systems usually 
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considered in this context, both the amplification and advection coefficients in (J6J) drop to zero as the material particle 
moves downstream through a stress field that smoothly approaches its free stream value. The slack pools behind the 
free stream, as shown in the lower inset of Figure [2] This should correspond to a location near the upstream base of 
the arch. 

OJ r 




FIG. 2. The grey curves are all ^ - 1, where a = ^- [l + tanh (^)] and r? = s + Tt. TOP AND CENTER: The evolution 
of a wave packet 8b(s, t) under (J3J (full black lines) and the integrated form of |5| (dashed black lines), for fi = T = 1, at times 
t — 10, 30, 50, 70, as viewed in a frame Eulerian on the string, with the downstream direction to the right and the horizontal 
axes representing s. The initial conditions at rj = s are 8b = 2bo for |3| and Sb — bo for jH]), where bo = ^sech (|^) cos(2ry). 

LOWER RIGHT INSET: The evolution of the slack variable / = \J (d s 8bf under (6} (full black lines) for fi = T = 1, at times 
t — 50, 150, 250, as viewed in a frame Eulerian on the string, with the downstream direction to the right and the horizontal 
axis representing r\ = s + Tt. The initial condition at r\ = s is uniform: I — . 

We have been talking about derivatives of the height rather than the height itself. However, in Figure [T] our 
slope amplifier is coupled to a height rectifier, namely a table. Is this sufficient for the formation of a structure? To 
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answer this, we have implemented a simple Verlet integration [55j of massive beads and stiff springs. This allows us 
to nullify bending and twist elasticity, frictional interactions with the table, and any other messy and experimentally 
unavoidable physics. The 5000 beads begin sinusoidally woven in a plane defined by the potential minimum created 
by a constant gravitational force perpendicular to a hard plane "table", an exponential with argument 100 at a 
penetration depth of one spring rest length. A velocity perpendicular to the weave is imposed on one end bead. A 
small quantity of uniform noise is added to the vertical force on each bead at all times, and a damping coefficient of 
a few percent is introduced during the first few timesteps to counteract errors in specifying initial bead positions. 

The snapshot of shape, stress, and height in Figure [3] is the result of a simulation with initial layout amplitude and 
frequency of ^6 cm and ^0.6 cycles/cm, parameters similar to Figure I] We have set /i = ||g|| = T = 1, and the 
half- width of the uniform noise distribution is 10 -3 . The extensions of the springs are used as a proxy for stress, and 
these data are spatially averaged. An initial attenuating pulse travels down the chain, leaving a nearly quiescent wake, 
followed by a rapid transition to the free-stream stress fiT 2 . This very noisy stress plateau also carries oscillations 
due to the undamped lateral motions of the chain arising from its initial configuration. An arch does indeed form, 
on the downstream side of the transition. The effect of gravity in the simulation is rather minimal, affecting only the 
qualitative shape of the arch. Without the "table" , however, the slack isn't rectified, and no arch forms. We also note 
that the dynamics are quite robust with respect to changes in the spring constant. The simulations shown in this 
letter involve strains of about 2.5% in the plateau, but the important features are quite similar for strains an order of 
magnitude larger or smaller. However, stiffcr springs result in proportionally larger oscillatory stresses with respect 
to the plateau stress. 
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FIG. 3. Height (grey line, normalized by the initial layout amplitude), stress (black dots), and averaged stress (black line, 
moving average over 50 beads, normalized by the theoretical free stream stress /iT 2 — 1) after 150 ms for a configuration with 
initial layout amplitude and frequency of ~6 cm and ~0.6 cycles/cm. The velocity boundary condition is imposed on bead 
#5000. A portion of the corresponding three-dimensional shape is shown in an oblique view in the upper left inset. Spring rest 
lengths are 2 mm, and the strain corresponding to a stress of 1 is 2.5%. The raw stress data extend from approximately -9 to 
+ 11 in the noisy plateau region. Noise from a uniform distribution of half-width 10 -3 is added to vertical forces; this may be 
compared to n\\g\\ = 1 and spring forces on the order of 10 7 . 



We return to our assumption of a steady stress distribution in the traveling wave frame. Figures [4] and [5] show 
snapshots of two chains at three identical times. The chain of Figure [4] is that of Figure |3j while that of Figure [5] 
corresponds to an initial arrangement with similar amplitude but lower frequency, ~6 cm and ~0.1 cycles/cm, and 
the same level of added noise. The high-frequency chain appears to sustain a steady, or perhaps very slowly growing, 
transition region between low and high velocity, and is highly receptive to noise. This region moves approximately 
6-7 beads/ms upstream, not much faster than the traveling wave coordinate (5 beads/ms). The transition region on 
the low-frequency chain steadily widens as the chain unfurls, its leading edge moving significantly faster upstream 
(9-11 beads/ms), and out-of-plane growth is slower to develop and less pronounced. Over the time interval shown, 
this region retains much qualitative information about its initial shape even as the material within it has achieved 



G 



highly tangential velocities. The stress on this chain appears to reflect this ordered arrangement, although the nature 
of the pulse train in Figure [5] remains to be explored. 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



stress , 
height 





1000 



2000 



3000 



4000 



bead tt 



5000 



5000 




bead tt 



5000 



FIG. 4. Height (grey lines, normalized by the initial layout amplitude) and averaged stress (black lines, moving average over 
50 beads, normalized by the theoretical free stream stress pT 2 = 1) after 55, 175, and 295 ms for a configuration with initial 
layout amplitude and frequency of ~6 cm and ~0.6 cycles/cm. Top views of a portion of the corresponding shapes are shown 
in the upper insets. Other details as for Figure [3] 
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FIG. 5. Height (grey lines, normalized by the initial layout amplitude) and averaged stress (black lines, moving average over 
50 beads, normalized by the theoretical free stream stress /xT 2 = 1) after 55, 175, and 295 ms for a configuration with initial 
layout amplitude and frequency of ~6 cm and ~0.1 cycles/cm. Top views of a portion of the corresponding shapes are shown 
in the upper insets. Other details as for Figure [3] 



In Figure [6j we show the maximum height of the high-frequency configuration of Figures [3] and |4j with levels of 
added noise between 10 -3 and 10 -15 . There is a clear trend towards an earlier appearance of the arch with higher 
noise levels, but even roundoff errors will generate nontrivial effects. 

Our preliminary study has raised several questions. Perhaps most interesting is the issue of how to analyze the 
planar dynamics we see in the insets of Figures [4] and [5j and the associated stress generation, in terms of the time- 
dependent screening and source terms of equation ([2| . There is some sort of signal being transmitted along our string, 
conveying information about the velocity boundary condition. The process of transmission is highly dependent on 
the evolving shape of the curve, and thus more complicated than transverse impulsive loading of a straight string 
[29) . Interesting physics arises from gradients in the stress, which are absent in the Klein-Gordon- type equations 
describing the perturbative dynamics around stationary states of elastic rods [30 , though we note that a uniform but 
time- varying tension is known to cause instabilities of planar motions of fixed-end strings [31] . What slowly evolving 
stress distributions can be sustained by moving elastic curves? And which spatially developing motions will inevitably 
bear noisy structures like our arch? Are there analogues in elastic sheets, nets, fluid membranes, or the Navier-Stokcs 
equations in certain geometries? Incompressibility of the system appears merely an analytical convenience, not a 
strict requirement. Finally, we have disregarded the periodic aspects of the curve's motion, and still lack an explicit 
description of the role of the rectifying potential. Is it possible to describe the arch growth as a secular result of some 
forcing by the table, perhaps on the high-curvature turning points of the curve, or is the resemblance to a resonance 
superficial? 
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FIG. 6. Maximum height (normalized by the initial layout amplitude) as a function of time for a configuration with initial 
layout amplitude and frequency of ~6 cm and ~0.6 cycles/cm, with noise added to vertical forces from uniform distributions 
of half- widths (from left to right or top to bottom) 10" 3 , 10" 6 , 10" 9 , 10~ 12 , and 10" 15 . Each curve is an average of three runs; 
take-off times of individual curves are highly reproducible. The lowest noise level is indistinguishable from roundoff noise up 
to a time of about 60-65 ms. Normalized spring rest lengths are ~0.033. 
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